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Abstract 

We apply recently developed integrable spin chain and dilatation operator techniques 
in order to compute the planar one-loop anomalous dimensions for certain operators 
containing a large number of scalar fields in M = 4 Super Yang-Mills. The first set of 
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operators, belonging to the SO(6) representations [J,L — 2 J, J], interpolate smoothly 
between the BMN case of two impurities (J = 2) and the extreme case where the 
number of impurities equals half the total number of fields (J = L/2). The result for 
this particular [J, 0, J] operator is smaller than the anomalous dimension derived by 
Frolov and Tseytlin |hep-t h/0304255 for a semiclassical string configuration which is 
the dual of a gauge invariant operator in the same representation. We then identify a 
second set of operators which also belong to [J, L — 2 J, J] representations, but which do 
not have a BMN limit. In this case the anomalous dimension of the [J, 0, J] operator 
does match the Frolov-Tseytlin prediction. We also show that the fluctuation spectra 
for this [J, 0, J] operator is consistent with the string prediction. 



*also at ITEP, Moscow, Russia. 



1 Introduction 



In the AdS 5 x S 5 conjecture [HI21E1, it was argued that the anomalous dimension of a 
generic gauge invariant operator scales as 
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(A) 1/4 v^, (1.1) 



where A = g^N is the 't Hooft coupling and i is the string level [2]. However, there 
are some operators which do not behave like (jl.ljl even for large values of A • O ne 
example are the half-BPS operators whose dimensions are protected by supersymmetry. 

Another example are the BMN operators, which are nearby to the half-BPS opera- 
tors 0. These operators have been heavily studied 000131^111111121113111111311211111 
IT8 l ll9 | l2n [ l2 H I22j. For these operators, one starts with a half-BPS operator Tr(Z J ) where 
% = 05+i06> and then inserts "impurities". The impurities, e.g. the scalars X = </>i + i</>2 
and y = 03 + z04, are inserted in various positions in the operator 

= TiZ...ZXZ...ZXZ...ZYZ...Z (1.2) 

Alternatively, one can think of this as an operator with L scalar fields, where L — J 
of the original Z's are converted to one of the other scalars. One-loop effects mix the 
operators among themselves and the operators with definite scaling dimension, along 
with their anomalous dimensions, are obtained by diagonalizing a matrix. This matrix 
of anomalous dimensions is particularly easy to compute by analyzing the action of the 
dilatation operator f3E]> even in the non-planar sector and at higher-loops. 

In the BMN limit, where the number of impurities is small compared to L, it was 
shown that the eigenstates of this dilatation operator are isomorphic to particles on a 
circle of length L, with the number of particles equal to the number of impurities. With 
a small number of impurities one can safely ignore the particle interactions and use the 
dilute gas approximation. The total momentum of the particles has to be zero in order 
to guarantee cyclicity of the trace, and one finds that the one-loop anomalous dimension 
behaves like 

7 = ^£^> = ^£w 2 (i-3) 

i i 

where pi are the individual particle momenta and e(pi) can be thought of as their energies. 
Quantization on a circle of length L then gives 

7=^5» 2 (1.4) 

i 

where the are integers. 

Remarkably, these results can be directly related to a string calculation on AdS 5 x S 5 

such that A' = X/L 2 is small. One then finds a one to one map of the BMN operators 
to string states propagating in the plane wave limit jSHEB] of AdS^ x S 5 . These states 
can be written as 
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where p + = L in appropriate units. In the limit of small A', the results then match with 
(Q. 

More recently, it has been suggested that other operators can be compared with 
semiclassical string calculations |39H40H4HI42|I43H44| . The operators of [4"2~ll4"4"j have order 
L impurities, with each impurity carrying momentum of order 1/L 2 . The anomalous 
dimension is now of order X/L. The corresponding string states are expected to have 
order L creation operators. 

In particular, Frolov and Tseytlin considered operators that have J' of the Z scalars 
and an equal number J of the X and Y scalars 44J. Assuming that these are highest 
weight states, one finds that these are operators in the [0, J' — J, 2 J] representation of 
SO(6) if J' > J, or the [J — J', 0, J + J'] representation if J > J', up to conjugation. In 
the extreme case where J' = 0, one finds a semiclassical result which is twice the BMN 
prediction That there is a difference from BMN is not surprising - the density of 
impurities is of order 1, so the dilute gas approximation is naively expected to break 
down. 

The case of a very large operator with a high density of impurities is not easily 
accessible by field theoretic means. Nevertheless we can still make progress in computing 
the one-loop planar anomalous dimension for these types of operators. This is because 
the planar one-loop dilatation operator for scalars, representing the interaction found 
in ^2j, is isomorphic to the Hamiltonian for an integrable system ^H] and one can then 
use all of the machinery that goes along with it. 

In particular, if one restricts the operators to be composed only of Z and X scalars, 
then the problem simplifies dramatically 1 . This is because the dilatation operator can 
be mapped directly to the Hamiltonian of the Heisenberg spin chain, which was solved 
long ago by Bethe. The Z scalars inside O are the up spins and the X scalars are the 
down spins. The field theory representation [J, L — 2 J, J] corresponds to a configuration 
of total spin L/2 — J. The eigenstates for the spin chain can be found by solving the 
Bethe equations. These are easy to solve for the case of two spin flips 2 , but for more 
than two flips these are in general only solvable numerically. 

Solving for a finite number of L/2 spin flips with total spin looks almost hopeless. 
However, in the large L limit, there is some hope of making progress since here one can 
convert the Bethe equations to an integral equation. This is how the ground state of 
the antiferromagnetic Heisenberg spin chain was solved; it has total spin and energy 
proportional to —L. 

To make contact with the proposal of 03] the low-lying states of the ferromagnetic 
spin chain with total spin and energy proportional to 1/L are of particular interest. 
Here, the antiferromagnetic ground state is the state with the highest energy. In contrast, 
the semiclassical string configuration is expected to correspond to the state of the lowest 
energy, or at least much lower energy. 

In this paper we solve this problem for [J, L — 2 J, J] representations where J and L 
are assumed to be large. This corresponds to operators of the form (modulo ordering of 

1 Operators with only Z and X scalars are special as they only mix amongst themselves due to charge 
conservation |21j . 

2 In the context of J\f = 4 SYM these are the finite J versions of the two-impurity BMN operators 
which were identified and diagonalized in [lfi| . 
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fields) 



O = Tr Z L - J X 



+ ... 



(1.6) 



As in the solution for the antiferromagnet ground state, this is solved by converting the 
Bethe equations to an integral equation. However, unlike the case of the antiferromag- 
netic ground state, the Bethe roots do not lie on the real line, but instead extend into 
the complex plane. 

We consider two classes of [J, L — 2 J, J] operators. The first class has two contours 
of Bethe roots symmetric about the imaginary axis which we call the double contour 
solution. The second class has all Bethe roots on the imaginary axis which we call 
the imaginary root solution. These latter Bethe states are "singular" in the sense that 
the Bethe equations need to be regularized. One of the main results of this paper is 
to identify the [J, 0, J] (spin 0) imaginary root solution as the SYM dual of the string 
solution found by Frolov and Tseytlin jUj. 

The agreement between the SYM and semiclassical string theory calculations is re- 
markable: the anomalous dimension exactly coincides with the energy of the string, 
and the spectrum of nearby operators coincides with the spectrum of small fluctuations 
around the classical string solution. Furthermore, the semiclassical analysis indicates 
that there is an instability in the string solution. This is consistent with our analy- 
sis, where we show that the imaginary root solution is not the lowest energy [J, 0, J] 
state. Instead we show that the double contour solution is the ground state for this 
representation 3 . 

We first consider the double contour Bethe state. In the limit of small J, this state 
approaches the BMN states with J/2 each of a_i and a + i oscillators. We are able to 
solve this problem by first considering the unphysical region where J < 0, and then 
analytically continuing to J > 0. The problem then reduces to solving for two elliptic 
equations which is easily done numerically to any desired precision. For the case where 
J = L/2, and in the strict J — > oo limit, we find that the anomalous dimension is given 



The result in (jl.7j) is bigger than the BMN result (jl.4j) where one finds (3 = 1/2. However, 
(jl.7|) is smaller than the semiclassical string result f3 = 1 in jH]. Hence, we conclude 
that this double contour solution is not the gauge dual to the semiclassical string. 

To help us find the gauge dual to the semiclassical string, we compare our result 
in (jl.7|) to actual numerical computations of [J, 0, J] ground states. We find very close 
agreement when J is even. However, when J is odd, we find that while the result in 
(jl.7|) appears correct, the configuration of Bethe roots is completely different from the 
even case. In this case there are two curving contours in the complex plane connected 
by a straight line of roots along the imaginary axis. It would be interesting to solve the 
integral equation describing this distribution of Bethe roots. 

3 Strictly speaking, the double contour solution is valid only for even J while the imaginary root 
solution is valid only for odd J. Indeed, the ground state for odd J has its Bethe roots in a completely 
different configuration than the double contour, and we have yet to determine the even J analog of the 
imaginary root solution. However, in the thermodynamic limit (large L), there should be no distinction 
between even and odd and we will provide strong evidence for this in Sec. [21 



by 




with /3 = 0.7120321458... 
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Nevertheless, it turns out that there are also solutions where all Bethe roots lie on the 
imaginary axis. With all imaginary roots, the relative coefficients in the sum over states 
for the spin chain (or gauge invariant operator) are real. For J small compared to L, the 
roots form two overlapping "strings" where the separation between roots is very close to 
i/2. If J is of the same order as L, then the outer roots spread out. If we then take the 
limit of large J and L, we can reduce the Bethe equations to an integral equation which 
is very similar to the equation found by Douglas and Kazakov for two-dimensional QCD 
on the sphere [3S1I3S]- We find that when J = L/2 there is a critical point and for this 
value the anomalous dimension matches the prediction in [33]. Unlike the case of the 
double contour, as the number of roots is reduced, the state does not approach a BMN 
state. 

To further verify that this imaginary root solution is the gauge dual of the semi- 
classical string, we consider spinless fluctuations about the solution. This is done by 
pulling roots off of the imaginary axis and onto (or close to) the real axis. Here we find 
a spectrum that is consistent with that found in |4*lj . 

In section two we solve for the solution with a double contour of roots in the large 
L limit. We find 7 as a function of L and show that it gives (|1.7)1 . In section three we 
compare this result to actual numerical computations of [J, 0, J] states. Here we show the 
surprising difference between the even and odd configurations. In section four we consider 
and solve for the imaginary root solution in the large L limit. We show that it matches 
the predictions in [33] for the anomalous dimension and the fluctuation spectrum. In 
section five we consider an analogous solution to an SO (6) singlet configuration and show 
that it matches a semiclassical prediction made in [32] • in section six we present our 
conclusions. 

2 The double contour solution 

In this section we describe the ground state for the [J, L — 2J, J] representations. Strictly 
speaking, this solution is valid only if J is even. If J is odd and L > 2 J, then the true 
ground state is doubly degenerate. But it approaches the double contour solution in the 
limit of large L. In the next section, we will argue that if J is odd and L = 2J, then 
the ground state will have a completely different root configuration than the solution 
described in this section. However, the anomalous dimension appears to be the same. 
We begin with the Bethe equations, which read 4 



Here L is the total length of the spin chain and J is the number of impurities 5 . The 
complex numbers Uj are the rapidities of the upside down spins (2uj = cot pj, where pj 
is the momentum of jth impurity). Normally the rapidities are real, but the impurities 

4 For a nice review see |47| . 

5 Note that this connotation has recently been introduced in the BMN literature. This is different 
from impurities in condensed matter theory, where one would rather speak, equivalently, of "excitations" , 
"particles" or "magnons". Or even simpler, "spins down"! 




(2.1) 
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can form bound states in which case their rapidities acquire imaginary parts. In the 
literature on the Bethe ansatz the bound states of impurities are usually referred to as 
"strings" . We shall study bound states of a macroscopically large number of impurities 
and shall argue that these large "strings" in the Heisenberg model are dual to strings in 
AdS§ x S 5 spinning around the five-sphere. The state with J impurities belongs to the 
representation [J, 2L — J, J] of SO (6) and we define the filling fraction 

J 1 , . 

a=- L <-. (2.2) 

We can also write down the logarithmic version of the Bethe equations 

Uj + i/2 -J-^ Uj — Uk + i / % 

L log -± L- = 2vrm, + V log k - : 2.3 

Uj—t/2 f— ^ Uj — Uk — % 

k±j 

where the numbers rij reflect the ambiguity of the branches of the logarithm. These 
numbers can be interpreted as the lattice momentum of the j-th 'particle'. 

Since the spin chain corresponds to an operator with a trace, we are only interested 
in states that are invariant under shifts. This forces us to set the trace condition 

n^f ■ 1 (2 ' 4) 

i=i J 

The eigenvalues of the planar dilatation operator can be found by solving the system of 
Bethe equations for the roots Uj (j = 1, . . . , J) [TH] . This is somewhat cumbersome for 
chains of small length, but very advantageous for long chains where a direct diagonal- 
ization [HH l2*T] is no longer feasible. The planar anomalous dimensions are then given to 
one-loop, for general L and J, by 

A = L + 7 with 7 =^£ lj (2.5) 

3=1 3 

In the "thermodynamic limit" where L — > oo, the roots are of order Uj = O(L). Let us 
therefore rescale them as 

Uj — ► L Uj (2.6) 
Then one easily proves that the Bethe equations ()2.3|) simplify for L — > oo to 

1 2 1 , . 

= 2-Krij + - > . 2.7 

Uj L ^ Uj - u k 

Here, the term on the left corresponds to a global potential, while the second term on 
the right corresponds to pairwise repulsion of roots. As is standard when taking the 
thermodynamic limit of the Bethe equations, eq. ()2.7)1 can be conveniently turned into 
an integral equation. 

In order to turn the Bethe equations ()2.7|) into a single integral equation we need 
to understand the structure of the distribution of Bethe roots in the complex plane. 
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Figure 1: Bethe roots. For large L the roots condense into two cuts 



To gain some intuition, let us concentrate on the case that L and J are large but the 
filling fraction a <C 1/2. We will also assume for now that J is even. The trickier 
case of odd J will be discussed (but not analyzed to the end) below. As can be seen 
from (J2.7jl . to a reasonable approximation the Bethe roots may be placed around the 
positions Ui = l/(27rnj), as long as no two n« are the same. However, we wish to lower 
the energy as much as possible, which means setting all rii to ±1 with an equal number 
of each to satisfy the condition in (|2.4jl . The interaction term in (|2.7|) pushes the roots 
apart. It is well known that the roots are pushed in the imaginary direction and form 
"strings" roughly parallel to the imaginary axis with the separation between the Ui of 
order 1/vZ [THj . It is clear from the Bethe equations that given a root «j, there has to 
be a root Uj = u*, therefore the distribution of roots is invariant under reflection about 
the real axis. 

It was also shown that there was a 1/L correction in the real direction towards the 
imaginary axis, where the correction is larger the further the root is from the real line. 
Finally, it was argued that placing the roots like this led to corrections to 7 which were 
of order 1/L 3 , while changing the branch gives a 1/L 2 correction. Hence, at least for 
small J/L, it is energetically favorable to evenly distribute the roots on the ±1 branches 
and hence the distribution looks like that in Fig. 

If J is of order L, then adding a root to the heavily populated ±1 branches costs an 
energy of order 1/L 2 . However, we will assume that for a ~ 1/2 it is still favorable to 
place the roots in these branches. We can also make a slight generalization and place all 
of the roots in the branches ±n. From the semiclassical string perspective, we expect 
this to correspond to a string that winds n times around itself. The net effect of this is 
to multiply vA, the effective string tension, by a factor of n. 

Assuming that J is still even, the roots are distributed so that they are symmetric 
under reflection about the imaginary axis. Hence we only need to consider the first J/2 
roots Uj with j = 1, . . . , J/2, since we can now assign Uj + j/ 2 = —Uj. Based on our 
previous argument, we now assume that rij = n = 1, 2, . . .. Let us introduce a density 



describing the distribution of roots with positive real part in the complex plane. This 
is expected to turn into a smooth function in the limit of large L, with a continuous 
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Figure 2: Cuts flip for a < 0. Both cuts can be mapped to one via the map u 2 = A + 



support tracing a curve C + (see Fig. Q). The density is normalized to one: 



du p{u) 



(2.9) 



c+ 



One easily verifies that the Bethe equations ()2.7|) can be reexpressed with the help of 
the density as 



dv p(v) 



1 

2a 



1 u 



(2.10) 



i c + u — v \za J a 

Here the integral has to be understood in the principal part sense, as we need to omit the 
case k = j in equation ()2.7|) . The anomalous dimension, cfeq. ()2.5|) becomes (remember 
the rescaling Uj — > Luj) 



A 



l = f~ 2 jl du^ (2.11) 
Sit 2 L J c + u 2 

Therefore all we need to do is invert the singular integral equation (|2.1()jl : this yields 
the density, and therefore, in light of eq. (|2.11jh the anomalous dimension. Luckily, 
similar equations have appeared previously in the case of large TV matrix models (which, 
incidentally, are also integrable models of sorts). In particular, equation 1)2.10)1 is a close 
relative of the equations describing the solution of Kostov's 0(n) multi- matrix model 6 . 
An inconvenient feature of ()2.1()jl is that the precise shape of the contour C + is unknown. 
The situation can be improved by first solving the equation in the unphysical regime 
a < 0, corresponding to a negative number of down spins. Here we expect the cuts to lie 
on the real axis, as is shown in the matrix model problems. The contour C + becomes the 
real interval [a, b] while the mirror contour C~ turns into the flipped interval [—b, —a], as 
in the left half of Fig. El This way we merely need to determine the endpoints a, b but 
not the shape of the contour. After solving the equation for general negative a we can 
subsequently analytically continue back to the physical regime a G [0, 1/2]. Next, let us 
apply the same folding map as in the case of the 0(±2) matrix model j48j : 

\{b 2 -a 2 ) 



p{u) 



u 



VA + Bx with A 



W + a 2 



This maps the two cuts [a, b] and [—6, 
[—1,1]. Our equation becomes, 

dyp{y)- 



and D = Ub--a-) (2.12) 
on top of each other, and to the interval 



x 



y 



- 2 



a. 



27m 



VA + Bx 



a 



(2.13) 



J see 



and references therein. 
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where p(u) = p(x(u)). Equation (j2.1Hj) is a standard one-matrix model (or airfoil) 
equation with a non-polynomial potential involving \/ A + Bx. The new cut C with 
branchpoints at x = oo and x = —A/B appears due to the above folding map, as shown 
in the right part of Fig. |21 The latter branchpoint lies for a < to the left of —1 on 
the negative real axis, while for a > it sits on the imaginary axis. The solution for the 
density is immediately obtained by an inverse finite Hilbert transform: 



p(x) 



1 VT 



vr 2 y/A + Bx J -i 
It simplifies to 



dy- 



27m 



y/A + By 



, , 2n yjl-x 2 
p[x) = — 



dy 



1 y/A + By 



(2.14) 



(2.15) 



na y/A + Bx J-i " x - y a/1 - y 2 

This is nearly the solution to our problem, except that we still need to determine the 
parameters A, B as a function of the filling fraction a. These are determined from the 
normalization condition ()2.9|) which reads in the new variables 



B 
~2 



p{x) 



and the further condition 



dx 



i dX \/A + Bx 



1 \ 2vm r- A — 

2 ] VA + Bx 

a J a 







VT^x~ 2 

which ensures the positivity of the density on the cut. It simplifies to 

1 



1 yfA- 

dx ■ 

l 



Bx 



VT 



x' 



— | - — a 

n 



(2.16) 



(2.17) 



(2.18) 



This completes the solution: A and B are determined through the two equations (|2.16|) . 
(|2.18jl . which one could express through elliptic integrals of the first and second kind, 
while the density is given by equation ()2.15j) . which could be found explicitly in terms 
of an elliptic integral of the third kind. However, as we are chiefly interested in the 
anomalous dimension ([2. II)) . we can simply eliminate the density. Plugging ()2.15|) into 
both the normalization condition ()2.16|) and the expression for the anomalous dimensions 
(|2.1ip . we respectively find, after exchanging orders of integration, 



and 



i 

dx —= 

i VT 

A 



x 



■VA + Bx 2nVA^ 



B^ 



7 



2a 



A 



(2.19) 



(2.20) 



16tt 2 L \y/A 2 -B 2 A 2 -B 2 / 

Thus 7 is explicitly given by the parameters A, B which in turn are implicit functions 
of a and n, as stated in eqs. (|2.18jl . (|2.19)1 . It is now a straightforward exercise to invert 
these equations in terms of a series in a to any desired order, one finds 
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Xn 2 J 
2L 2 



(1 + \a + |« 2 + §« 3 + ifa 4 + + §t^-- •) (2-21) 



S 



While we derived this result for negative a, we see that the function is perfectly analytic 
at a = and we can immediately analytically continue to the correct Bethe phase a G 
[0, 1/2]. The leading order 0(a) result correctly reproduces the dilute gas approximation 
(II. 4J) and the famous BMN result Xn 2 /L 2 for J = 2. Let us also note that the leading 
1/L correction to the BMN limit 

was obtained in it agrees with the 0(a 2 ) term in (j2.21|) and thus yields a nice 
check of our method. The above series ()2.21|) converges rather rapidly in the interval 
a G [0, 1/2] allowing us to obtain the anomalous dimension to arbitrary precision for any 
a. We can also expand the anomalous dimension around a = 1/2, where we find 

7 = [0.7120 - 1.0916(1 - 2a) + . . .] (2.23) 

A plot of these results in presented in Fig. |3] in the next section. 

An alternative method for solving ()2.18|) and ()2.19|) is to analytically continue a first, 
thus moving B onto the imaginary axis. For small values of a, A is positive, but decreases 
as a is increased and eventually changes sign for a critical value of a. At this value the 
branch point on the imaginary axis touches the cut along the real line. 

To avoid this collision, we simply deform the integration on the real line to an inte- 
gration along the cut on the imaginary axis. If we define C = —iB, then the equation 
in ()2.19|) becomes 

2 / dy ■ ■ = ; 2.24 

Ja/c ^TT^VCy - A 2n VA 2 + C 2 

Defining a new variable uj = \JCy — A, ()2.24|) can be reexpressed as 



oo 1 1 1 

du = ; . (2.25) 

-oo 



v/(cj 2 + A) 2 + C 2 2n VA 2 + C 2 



We then do the same with (|2.18|h although now there will be a divergent piece that needs 
to be subtracted off. Deforming the contour, redefining variables and subtracting off the 
infinite piece, we find 

°°do;(l--= u2 )=i(l-g). (2.26) 

oo \ v^ 2 + A) 2 + C 2 J n\2 J 

If we set A to zero, then the elliptic integrals reduce to ordinary integrals and we can 
find an analytic solution. In this case 

a = - - - 4?1 " x ^ w 0.2715 < -, 
2 (r(l/4)) 4 2' 

\n 2 

^si- (2 ' 27) 
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Hence A becomes negative before a reaches a = 1/2. Assuming that A < and letting 
7] = —C/A we define the integrals 




(2.28) 



We then find from (l2~231) and (l2~2T)J) that 



a=;|l--^£^). (2.29) 

While it is possible to solve this for all physical values of a, let us only consider the case 
a = 1/2, which corresponds to setting I\ to zero. This is easily done numerically, where 
one finds 77 = 1.16220056. From this it is straightforward to show that the value for 7 is 
that in (jl.7j) and that the values for A and C are 

A = -0.086987288 C = 0.10109667. (2.30) 

In terms of our original variable u, we see then that the branch points are at 

u = ±0.1523 ±0.3319 i. (2.31) 

The exact distribution of the roots is found by starting from the branch points and 
finding the contour where du p{u) is positive definite. The positivity of the density 
unambiguously determines the shape of the contour jlHj. Here, p{u) is the analytic 
continuation of (|2.15jl . In terms of x, this means 

BP{X) dx (2.32) 



y/A + Bx 



is positive definite. While we have not solved for this over the whole contour, we have 
determined that the slope of the contour away from the upper right branch point in 
(|2.31j) is —0.743. Hence the distribution has the general form in Fig. ^ 



3 Comparison of results 

We now compare our results of the previous section to the exact anomalous dimensions 
of operators of the form (jl.6|) for reasonably large values of L. These can be obtained 
directly using the dilatation operator or using the exact Bethe equations. 

3.1 Dilatation operator 

For given L and J we collect all operators of the form Tr Z L ~ J X J ± . . . and act on them 
with the dilatation operator |21j . In this way the matrix of anomalous dimensions is 
generated and we subsequently diagonalize it. A slight complication is that the set of 
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J 7 [AJ/L 2 ] 


J 7 [AJ/L 2 ] 


2 0.6079271019 
4 0.6576949106 
6 0.6759273535 
8 0.6848144787 


3 0.9118906528 
5 0.8416441952 
7 0.8155060662 
9 0.7973649389 



Table 1: Numerical values for the lowest eigenvalues in [J, 0, J] 



operators (jl.6|) comprises not only the operators in the representation [J, L — 2 J, J], but 
also the ones for smaller values of J. These additional eigenvalues have to be removed as 
some of them are generically smaller than the one we are interested in. An easy way to 
achieve this is to repeat the analysis for the operators Tr Z L ~ J+1 X + . . . and remove 
the corresponding eigenvalues 7 . What remains are exactly the anomalous dimensions of 
the operators in the representation [J, L — 2 J, J] 8 . 

In this way we have obtained the anomalous dimensions for all operators of length 
L < 19. We give the numerical values for the lowest eigenvalue operators in Tab. ^and 
Fig. El We note that the behavior of the anomalous dimensions is significantly different 
for even and odd values of J. We will therefore treat them separately in what follows. 

Already for J = 8 the accuracy is very good compared to the numerically exact result 
()1.7|) at J — > oo. We can improve the estimate for J — > oo by linearly expanding around 
J = oo, 

T =(a + 6/J+...)^ (3.1) 
We fit the coefficients a, b to the available data for J = 4, 6, 8 and get 

7= (0.712-0.22/7+...)^ (3.2) 

The estimates depend on the selected set of data and the fit function, overall the agree- 
ment with ()1.7|) seems to be within 1%. 

Unlike the even case, we see from Tab. ^that for odd values of J the energies seem 
to approach an asymptotic value of O^L^/J) from above. A fit for odd J with (13. lj) 
yields 7 = (0.743 + 0.49/ J) A J/L 2 . While the even and odd case seems to asymptote 
to a different value, this could be due to the paucity of data points in the fit. In the 
thermodynamic limit where J is large, it is hard to believe there should be a distinction 
between even and odd. Further down in this section, we will present evidence that indeed 
the difference vanishes, see ()3.9|) . 

We can make further observations on the obtained data for L < 19. An important 
feature of the one-loop planar dilatation generator or of integrable spin chains in general 
is that states may come in pairs, see |21j . We observe that 

• the energy of the ground state is 0(1/L), 

• for even J the ground state is unpaired, 

7 One of the benefits of the Bethe ansatz is that states of a lower representation are easily identified, 
they have roots at infinity and discarded. 

8 For the representation [9,1,9] of SU(4) this involved 4862 states. Out of these, 884 belong to the 
representation [9, 1, 9], the other belong to [J, 19 — 2J, J] with J < 9. 
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Figure 3: States of lowest energy in [J, L — 2 J, J], J even. The plot shows how the energy per 
spin flip increases from the dilute gas, a = 0, to the maximum filling a = 1/2. Discrete values 
for L < 19 are obtained as eigenvalues of the dilatation generator and the curve for L — > oo 
represents the solution to the Bethe equations. The states with J = 2 were found in [TB] . 

• for odd J the ground state is paired unless J = L/2, 

• for odd J the energy of the lowest unpaired state is 0(1/ J), 

• for odd J and L all states are paired. 

It is also interesting to investigate the operators of lowest energy. For a = 1/2 the 
operators are 

j-i 

O = aTrZ J X J + b >n,n Tr Z m X n Z J - m X J ~ n + . . . (3.3) 

m,n=l 

where a and b are mixing coefficients. We find that a is by far the most dominant 
coefficient. The coefficients b m ^ n are less important, while the remaining coefficients are 
almost negligible. This is in analogy to the Weiss domains configuration of a ferromagnet. 
We find that alike spins cluster in domains with only the domain walls contributing to 
the energy. For a minimal energy solution the number of domains must therefore be as 
small as possible. 



3.2 Bethe roots 

The Bethe equations are non-linear equations of multiple variables Uj. It is therefore a 
difficult task to find solutions to them. A numerical solution requires the knowledge of 
the approximate position of the roots Uj. It would therefore be useful to know how the 
roots arrange in general. We have worked out solutions for small values of J and L to 
obtain some intuition in this respect. 

First of all we note that the Bethe equations (|2.1|) . the trace condition (|2.4jl and the 
energy ()2.5|) are invariant under negation of roots Uj — > —Uj. Therefore each solution 
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{uj} has a partner {— Uj} of equal energy unless {uj} = {— Uj}. This is the same pairing 
discussed in the previous subsection. Here we will concentrate on the unpaired states. 
This is also simplifies the analysis because now there are only half as many equations to 
solve. 

For unpaired solutions we must distinguish between an even and an odd number of 
roots J. For even J all roots come in pairs uj/ 2 +j = —Uj and the trace condition ()2.4|) is 
automatically satisfied. The important feature for odd J is that one of the roots must be 
zero, while the others come in pairs. At first glance, this would seem to violate the trace 
condition since the zero Bethe root contributes a factor —1 while Uj and — Uj together 
contribute a factor of 1. 

There is one way around this problem, which is to also place two Bethe roots at the 
singular positions ±z/2. These roots need to be regularized since the Bethe equations 
diverge for these values 9 . Hence we may assume that there are three roots 

uj = 0, mj-i.j-2 = i(±l/2 +e±6), (3.4) 

where 5 = 0(e L ). The Bethe equations associated to the singular roots and the trace 
condition are satisfied in the limit e — > 0. The Bethe equation associated to the zero 
root implies that L must be even. This is in agreement with the observation of the last 
subsection. The contribution of the special roots to 7 is finite and is given by 

A , / 1 1 \ 3A , . 

7 = — lim 4 - I = - — - . 3.5 

8ir 2 £^o\ e + e 2 -e + e 2 J An 2 y ' 

If J = 3 then there are no other roots and this state corresponds to the special three 
impurity operators discussed in fH\ . 

We have determined the root configurations for all unpaired states of the [J, 0, J] 
spin chain for 2 < J < 7. The positions of the roots were found by rewriting the Bethe 
equations in polynomial form. These equations were combined into a single polynomial 
equation of one variable using the resultant. All Bethe roots are then among the solutions 
to this equation. We summarize our findings in Tab. Eland Tab. El El at the end of this 
paper. 

We make the following observations. All the roots of a given mode number n form a 
vertical string of roots approximately at the real coordinate L/27in. This is in agreement 
with ^H]. What is different is the distance of roots in the imaginary direction. For a low 
density of roots a, the roots of equal mode number n should be separated by iO{yL). 
For a high density the separation is very close to i if the roots are close to the real axis, 
but spreads out as we move away from the axis. 

3.3 The odd ground state 

Since the three roots in ()3.4j) contribute a finite amount to 7, the other roots must 
somehow cancel off most of this contribution to leave a 7 of order 0(1/ L). We note that 
if a Bethe root is on the imaginary axis with \ui\ > 1/2, then this root contributes a 
negative amount to 7. 

9 The existence of these singular solutions to the Bethe equations was recognized by Bethe himself [50| 
and has drawn certain attention with regard to the completeness of the Bethe ansatz pi ,52, 53 ,53]- 
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J 


7 [AJ/L 2 ] 


Bethe roots 


n 


2 


0.607927 


0.288675 


1 


3 


0.911891 


±i/2 


2* 


4 


0.657695 


(±0.463265 ± 0.502294i) 




1 1 




1.104851 


±1.025705* ±0.041309 




03 




2.290302 


±0.525012 ±0.129473 




13 


5 


0.841644 


±i/2 


±0.998506i 


3*1 




1.291310 


±i/2 


±1.570673i 


4*0 




2.289413 


±i/2 


±0.638965 


4*1 




3.176721 


±i/2 


±0.236124 


4*2 



Table 2: States of [J, 0, J] 



Let us examine more closely the Bethe roots for relatively small odd J. Assuming 
the presence of the singular roots in (|3.4p . the [5,0,5] state has two more roots which 
are the negative of each other and which have to be either real or imaginary in order to 
be invariant under the shift operator. Hence, the Bethe equations in ()2.1|) reduce to a 
single equation 

( u + i/2 \ 8 _ (u + i)(u + 3i/2) 

\u-i/2J (u-i)(u-3i/2)' 1 J 

One then finds that the solution that gives the entry in Tab. ^ is imaginary and very 
close to i. A similar analysis for the ground state of the [7, 0, 7] operator shows that two 
more roots are placed very close to ±3z/2 and that the other two roots move even closer 
to ±i. This pattern continues up to and including the [11,0,11] operator, where roots 
are added near half integer imaginary values and where the lower roots move closer to 
the half integers as higher roots are added. 

The half-integer periodicity of imaginary roots can be explained by the same argument 
that leads to the string hypothesis for the structure of Bethe states in the thermodynamic 
limit |47j. If u n = iq n , then the Bethe equations in ()2.1|) take the form 




For q n of order one and positive, the exponentially large factor on the l.h.s. of ()3.7|) has 
to be compensated by a small denominator on the r.h.s., which arises when two roots 
are separated by an amount which is close to 1. Since we are starting with the roots in 
()3.4|) . this tends to put the new roots close to the half-integers. We can also see that 
increasing L forces the lower roots closer to the half-integers. 

One can also see that with this configuration of roots 7 will be of order 0(1/ L). If we 
assume that the roots are equally distributed along half-integer imaginary values, then 
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Uk ~ ki/2 

Rc u 



Figure 4: Distribution of roots for the odd ground state. 



from ()2.5|) we find 





J 1 [\J/L^\ 


J 7 [AJ/L 2 ] 


3 0.911891 
5 0.841644 
7 0.815506 
9 0.797365 
11 0.783974 


13 0.773992 
15 0.766364 
17 0.760374 
19 0.755558 
21 0.751605 


23 0.748305 



Table 3: Odd ground state energies for [J, 0, J] 



7 



8tt 2 



6-8 



J=2 



n- 



8vr 2 
8 A 



6-8 



1 1 \ 

3=2 j= i±l I 



(3.8) 



The prefactor in ()3.8J) is quite close to the values in Tab. It is also in agreement with 
the observation of Sec. 13. II that the energy of odd unpaired solutions is of order 0(1/ J). 

However, for large enough J and q n we expect this picture to break down since the 
l.h.s. of ()3.7j) will then be of order 1. In fact, for the [13, 0, 13] case we found that the 
ground state has 9 roots on the imaginary axis and 4 split off from it into the complex 
plane. For higher values of J, more roots will break off from the imaginary axis, see 
Fig. HJ In Tab. Elwe show the ground state values for 7 up to J = 23. If we fit the higher 
data points starting with J = 13 to the curve in (|3.1j) we find 



7= (0.715 + 0.77/ J) ^ 



(3.9) 



indicating that the odd result is asymptoting to (jl.7j) . Hence the even and odd cases 
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J 


7V l^-J/L 2 } 


72,j [AJ/L 2 ] 


7s,j [AJ/L 2 ] 


3 


0.911891 






5 


0.841644 


1.291309 




7 


0.815506 


1.033470 


1.684344 


9 


0.797365 


0.963623 


1.210401 


11 


0.783974 


0.931003 


1.068883 


13 




0.907919 


1.008200 


15 




0.888306 


0.976885 



Table 4: [J, 0, J] states with roots on the imaginary axis 

are approaching the same result even though the configuration of Bethe roots is totally 
different. 

3.4 States with imaginary roots 

Even though they will not give the smallest value for 7, there still can be solutions to the 
Bethe equations with all roots on the imaginary axis. These remain interesting objects, 
because they still correspond to some string states in AdS§ x S b . What one finds is that 
the lower roots are very close to the half-integers, but the higher roots tend to spread 
out. In fact, for a given value of J > 3, there seem to be multiple solutions of imaginary 
Bethe roots. Results for 7 where the state has all imaginary roots are shown in Tab. 0] 

Let us consider all imaginary solutions for a given value of J, we label them Uk,j- 
Comparing this to the case J + 2, it seems that for each solution Uk y j there is a cor- 
responding solution Uk t j+2- I n the examples considered in Tab. 0] this solution has a 
similar distribution of roots with two additional roots. Often, these roots are also on the 
imaginary axis, in some cases, however, they can spread out into the complex plane. For 
example, this happens in Ui^, see the discussion in the last subsection. Although the 
roots cease to be exclusively on the imaginary axis, the energy 7113 continues to follow 
the extrapolated trajectory, see Tab. El However, when we discard these states we notice 
that the value of 7L 2 / J jumps up, see Tab. HI 

We can therefore assume the following picture: Trajectories of imaginary solutions 
Uk,j appear at some value of J. At a higher value of J some roots of Uk,j split off of 
the imaginary axis. When we discard such states, which have a relatively low energy, 
the state with the minimum energy will be Uk+i,j and the lowest jL 2 / J lurches upward. 
But one can also see in Tab. 0] that as J becomes large, there is some indication that 
the result is approaching the semiclassical result of Frolov and Tseytlin for the spinning 
string. We show that this is the case in the following section. 

4 The gauge dual of the semiclassical string 

In this section, we show that the state with all imaginary roots has the expected prop- 
erties of the dual to the semiclassical string in We first show that the anomalous 
dimension has the expected behavior in the limit of large L. We then consider spinless 
fluctuations about the solution and show that they too have the expected behavior. 
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4.1 A solution for imaginary roots in the large L limit 

In the previous section, we saw evidence that if all roots have the form Uj = iqjL with 
qj real, then the roots close to the real line are very close to the imaginary half integers, 
but those further away begin to spread out. We will call those roots separated by half 
integers the "condensate" . When L becomes large ()2.Hj) becomes 

M E — <«> 

where (14.1)1 is valid for qj outside the condensed region. Note in particular that there 
is no integer on the r.h.s. due to a branch ambiguity. The equation in ()4.1j) can be 
rewritten as an integral equation 

i=/^, (4.2) 
q J q-q' 

which is valid for q outside the condensate and where the root density cr(q') is given by 

*(90 = f E stf-n)- ( 4 - 3 ) 

k=-^l 

We now make the ansatz that the roots are condensed in the interior region, but 
start spreading out at some point of order L from the origin. This means that the root 
density is of the form 

— s < q < s, 
q < —t or t < q 

where t and s are to be determined. Notice that this configuration is very similar to 
one considered by Douglas and Kazakov for the strongly coupled phase of QCD 2 on the 
sphere (25!) where the density also has a condensed region. The problem considered 
here is somewhat different since the QCD 2 equations would have q on the l.h.s. of (|4.2jl 
instead of \jq. In terms of a(q), (|4.2|) becomes 

I_41 n i±£= fdq'aiq') (—!— + —!—) . (4.5) 

q q-s j s \q-q' q + q'J 

The solution of this integral equation can be written in the form 

a(q) = -V(? 2 - s 2 )(t 2 - q 2 ) I -— + 4 [' . = =] . (4.6) 

If we plug (|4.fjjl back into (|4.5|l . we see that s and t must satisfy a consistency condition 

dV 1 (4.7) 




v/(s 2 -v 2 )(t 2 -v 2 ) st 
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We also have the normalization condition for a(q) 

j dq a(q) = 8s + 2 J dq a(q) = 2a , (4.8) 

where a is the previously defined filling factor. If we now insert (|4.fi|) into ()4.8|) and use 
(14. 7 j) we find the equation 

dv v 2 , . A . 

1 - 2a . (4.9) 



y/(s 2 -v 2 )(t 2 -v 2 ) 

The anomalous dimension is found by using (|2.5jl . ()3.8|) and (|4.4|) . which gives 
A 



7 



8tt 



6 - 2 S^-l?i) +0(1/L2 



" ^-if*Sj +0 (i/ ir (4.10) 



8tt 2 L 7 s v 2 

The integral inside ()4.10|) is found by using (|4.6|) . ()4.7|) and (j4.9j) . which leads to 

A /l 1 l-2a\ 

Let us now focus on the limiting case a = 1/2 which corresponds to the [J, 0, J] 
representation. Using ()4.9j) and (|4.7|) . we see that in the limit a — > 1/2, t and s approach 

t -> oo s -> — . (4.12) 

4vr v y 

Therefore, setting a = 1/2, we find 

which is precisely the string theory prediction in |44j . Starting from this result we can 
also expand the energy around a = 1/2. We get 

7 = A (1 + g(i - a) 2 + 24(| - a) 4 + 96(| - a) 6 + 408(| - a) 8 + . . .) (4.14) 

Interestingly, the energy is symmetric around a — 1/2 and reaches a maximum at a = 
1/2. 

We can also see that the solutions for the [J, L—2J, J] representations do not approach 
a BMN limit when J ^ L. The easiest way to see this is to realize that the roots are 
all close to imaginary half- integers and so we can borrow the result in (|H.8|) . giving 10 
7 = 2X/ir 2 J if J L. Hence, for finite J the anomalous dimension is finite in the large 
L limit. 



10 The same results can be derived from eqs. Ij4.7l) . 1)4. 8[l . (|4.11)l by noticing that s m t w a/4 and 
i — s ~ exp(— 4/a 2 ) as a — > 0. 
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Note also that for a = 1/2, the root density simplifies to 



?(,) = 4 ( 1_ /~5) ' (415) 

Hence, 4 — cr(q) has a Wigner distribution in 1/g. This is like the critical point for 
Douglas-Kazakov, where the eigenvalue distribution changes to a Wigner distribution at 
the critical area. 



4.2 Fluctuations 

We can also show that the fluctuations of the solution in the previous section are consis- 
tent with the Frolov-Tseytlin prediction. The fluctuations we consider are spinless, that 
is, they do not change the SO (6) representation. Hence, these fluctuations are found 
by moving Bethe roots around in the complex plane, but not actually changing the net 
total of roots. 

With this in mind, let us suppose that two roots are moved from the imaginary axis. 
If we assume that the remaining roots stay on the axis then in order to satisfy the trace 
condition, the two roots must be at ±/i where \x is positive real. In the large L limit, /i 
satisfies the equation 

1 , f ,„ <*) 



2im + / dq v ^ , (4.16) 
IX J n-iq 

where we have again rescaled the roots by a factor of L. The integer n corresponds to 
the different branch choice for the log. The root at —fi is on the — n branch. We may 
also assume that the function a(q) is the same one computed in the previous section. 
Using a(q) = cr(—q) and the expression in ()4.15jh we obtain the equation 
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, 2m+ L dq ^- il \L iq+ L dq )Z^rjk (417) 

Deforming the contour, we then find 

fi- 1 = 27iVn(n + 4) . (4.18) 

Hence, we see that n > for fi to be positive real. The contribution of these roots to 
the anomalous dimension is 

7m = 8tt 2 LV = ( L 2 } ■ (4 ' 19) 

The roots at ±/i also back react on the imaginary axis roots. To find this effect, we 
note that the roots at fi modify the equation in (|4.5j) to 

l-41ni±^- ' 9 2 = S t dq ' Hq > ) (l+^-) . (4.20) 
q q-s L(g2 + yU 2) J g {q } \q - q< q + q<) 1 } 

We can then solve for a(q) with the new term on the l.h.s. of the equation in ()4.20|) . 
Thus, we have 

r* 4 

dq a(q) = 8s + 2 dq a{q) = 2a - -, (4.21) 
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s 



dv 



1 4 1 



4 



y/( S 2 -V 2 )(t 2 



V 2 ) St L ^( S 2 +/i 2)( t 2 +/i 2) ' 



(4.22) 



A' 



and the same equation in ()4.9j) . Hence, we see that the distribution still has t — > oo if 
a = 1/2, but the position of s is shifted to 



We now compute the contribution of the imaginary roots to the anomalous dimension, 
where we find 



If we now let a = 1/2, use (ETTKj) and (pT23j) . and then add (ETTH1) to (ET241 . we find that 
the change in the anomalous dimension from ()4.13|) is 



This is the prediction of Frolov and Tseytlin, where our definition of n is shifted by 2 from 
theirs, and with extra factor of 2 since this configuration is essentially two fluctuation 
quanta. 

In their semiclassical analysis, Frolov and Tseytlin also identified an unstable mode, 
that is, a mode with an imaginary frequency. Since we are considering eigenstates of a 
Hermitian Hamiltonian we will not see this directly. It would seem that A7 is imaginary 
if n = —1 or —3. However, in this case fi would also be imaginary. But the derivation 
of ()4.18|) assumed that Re(/^) > 0. Hence the imaginary /i solution is not actually a 
solution to the Bethe equations. 

On the other hand, we have already established that there are states with lower 
energy and in the large L limit we would expect them to form a continuum, see Sec. 13.41 
Hence there can be a superposition of states centered about the imaginary root solution 
that will decay. If the superposition is spread out over solutions with width 0(1/L 2 ), 
then this will lead to frequencies with imaginary parts of the same order. 

Even though the imaginary A7 do not correspond to solutions to the Bethe equations, 
they still give the predicted imaginary frequencies of Frolov and Tseytlin. Hence a 
possible interpretation is that the states constructed with these values of /1 are precisely 
these superpositions of energy eigenstates. 

Placing roots at ±// is consistent with level matching of the fluctuations. If we were to 
consider a finite number of roots removed from the axis and placed at ±/^, we may ignore 
their interactions among themselves and so the change in the anomalous dimension is 
additive. If we were to put two roots on the same cut n, they would split off from the 
real line [TBI, but the splitting is of order 1/yL and so this can also be ignored at order 




(4.23) 




(4.24) 



A 7 = — (n + 2)^n(n + 4) . 



(4.25) 



1/L 2 . 
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From the string perspective, there can be nonsymmetric configurations of roots that 
are still consistent with level matching. Assuming that such states exist, then the imagi- 
nary roots will also be nonsymmetric, which means that they must stray from the imag- 
inary axis. It would be interesting to find these solutions. Also it would be interesting 
to find an analog of the semiclassical string configuration for even J. 



5 An SO (6) singlet solution 

It turns out that we can find a solution for an SO (6) singlet that is very similar to the 
solution in the last section. Once we consider the full SO (6), then it is necessary to 
introduce two other types of Bethe roots 18 . If we assume that all roots are imaginary 
and that the distribution of the two new sets of roots are the same, then the full SO (6) 
Bethe equations become 

q n + l/2\ L I r Qn Qm 1 i r ( Qn Tk 1/2 



TT Hn ~ Hm T 1 TT 

q n -l/2) 11 q n -q m -l 11 \q n -r k + 1/2 



n r J ~ r k + 1 TT r_j -q m - 1/2 
, ( .r i -r fc -lllr i - 9m + l/2 1 ' ' 

k^j J m J ' 



where rj refer to the new set of roots. Taking logs, rescaling by L and converting these 
equations to integral equations, we find 



q J q — v J q — v 

0= Idv^l-Udv^l (5.2) 
/ r — v 2/ r — v 



where a(q) is the same as ()4.3|) and u(q) is 



2 x - 

w(r) = — 2_^S(r - r k ) . (5.3) 



2 r d v^ (5.4) 



Therefore, we have 

q J q — v 

and so we have almost the same equation as in the previous section. The factor of 2 
leads to the modified equations 



4 



dv 2 



^/(s 2 -v 2 )(t 2 -v 2 ) st 
dv v 2 



-v 2 )(t 2 -v 2 



2 - 2a 



^ = T^{7 2 + T 2 - 2 —) • (5 - 5) 

In the SO (6) integrable chain there can be as many as L of the q type roots and 
L/2 of the r type. In this case, the representation is a singlet. Hence we see that the 
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maximum filling fraction is a = 1, which is reflected in the r.h.s. of (|5.5jl . We also see 
that s = l/(2ir) when a = 1. Therefore, the anomalous dimension is 

1= IZ ' (5 ' 6) 



which matches the semiclassical prediction for a circular string pulsating on 5*5 [42 . 



6 Conclusions 

In this paper we have constructed two types of gauge invariant operators which are in 
[J, L — 2 J, J] representations. The first type is a limit of BMN operators and was shown 
to have the smallest anomalous dimension for this representation. The second type of 
solution had a higher anomalous dimension, but was shown to have the same anomalous 
dimension and fluctuation spectrum as the semiclassical string, strongly suggesting that 
this is the gauge dual. 

The fact that the semiclassical string result in ^1] does not correspond to the operator 
with smallest anomalous dimension is mildly surprising; one usually expects a state with 
a large amount of symmetry to have a low, if not the lowest energy. It suggests that there 
are other semiclassical solutions to be found. It also suggests that these semiclassical 
solutions will involve elliptic integrals in order to reproduce the anomalous dimension. 

It would be very important, in order to check and complete our ideas, to find the 
solution of the Bethe integral equation for the ground state in the case of odd J, c.f. sec- 
tion IrS.rSl Likewise, it would be exciting to find the analog of the purely imaginary 
solutions for even J. On physical grounds, we expect no difference in the energies be- 
tween even and odd for large J, but as we have seen, the root distributions are certainly 
very different. 

The computations presented here are one-loop calculations, but nevertheless not ones 
that would be accessible to standard Feynman diagram and combinatorial techniques, as 
the real difficulty is the diagonalization of the states. This nicely illustrates the power of 
using the Bethe ansatz to do one-loop computations of anomalous dimensions in SYM 
and suggests that it might be applicable to even more gauge invariant operators. It also 
illustrates the extremely rich mathematical structure hidden in M = 4. 

Finally, let us stress the fact that the problems solved in this paper are quite different 
from the ones considered by condensed matter theorists. For magnetic chains one is 
either interested in the ferromagnetic or the antiferromagnetic phase. The Bethe vacuum 
(corresponding to the BPS states in gauge theory) is the true vacuum of the ferromagnet, 
but only a "fake" vacuum of the antiferromagnet. For finite filling fraction a, the two- 
contour states we considered are hybrids - they are "antiferromagnetic" in the sense that 
they have low spin, but ferromagnetic in the sense that they have the smallest energy, 
i.e. anomalous dimension, for fixed a. For this reason and to the best of our knowledge, 
despite all the work on the XXX Heisenberg chain since its invention more than 70 years 
ago, this problem has not previously been solved. 

Note Added: After this paper was completed we learned that Frolov and Tseytlin 
have studied fluctuations around the stable counterparts of the string soliton dual to 
the Bethe states discussed here. They found that the 0(X/L) term in the anomalous 
dimension is indeed not renormalized by the string corrections [53] . 
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Second Note Added: After this paper was submitted to the arXiv we were informed 
that Frolov and Tseytlin were able to find classical string configurations which are dual 
to the double contour solution presented here. |BT)| . 
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Tables 



-f[\J/L 2 ] 


Bethe roots 


n 


0.675927 


(±0.676245 ± 0.993633?) 


±0.678017 


111 


0.948774 


(±0.120258 ± 0.500000?) 


±1.638344? 


220 


1.485952 


±2.177728? ±1.000078? 


±0.001879 


005 


1.913430 


(±0.347918 ± 0.5000192) 


±0.695172 


221 


Z.oUoZOU 


±1.003317?' ±0.805340 


±0.012646 


U 1 


2.804933 


(±0.519954 ±0.500990?') 


±0.279477 


114 


3.163442 


±1.012980? ±0.386936 


±0.024593 


025 


3.283924 


(±0.560191 ±0.501721?) 


±0.083335 


115 


3.854384 


±1.063338? ±0.186496 


±0.047508 


035 


5.098922 


±0.657299 ±0.282367 


±0.084694 


135 



Table 5: Unpaired states of [6, 0, 6] 



7 [AJ/L 2 ] 


Bethe roots 


n 





815506 





±?/2 


±1. 000016?' 


±1.478037?' 


6*11 


1 


033467 





±?/2 


±0.999996? 


±2.241154?' 


5*10 


1 


684344 





±?/2 


±1.500899? 


±2.862938? 


6*00 


1 


870623 





±?/2 


±0.999885?' 


±0.882971 


5*11 


2 


366246 





±?/2 


±1.524153?' 


±0.985655 


6*01 


2 


749876 





±?/2 


±0.999612?' 


±0.449560 


5*12 


3 


065353 





±?/2 


(±0.701060 


± 0.504190?') 


6*11 


3 


156982 





±?/2 


±1.555226? 


±0.522436 


6*02 


3 


514967 





±?/2 


±0.998615?' 


±0.231443 


5*13 


3 


915465 





±?/2 


±1.584609? 


±0.295822 


6*03 


3 


939509 





±?/2 


±0.986006? 


±0.062156 


5*14 


4 


462322 





±?/2 


±1.603217?' 


±0.136490 


6*04 


4 


926542 





±?/2 


±0.740432 


±0.348252 


6*13 


5 


587007 





±?/2 


±0.752202 


±0.154446 


6*14 


6 


559173 





±?/2 


±0.364262 


±0.158450 


6*24 



Table 6: Unpaired states of [7, 0, 7] 
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